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Abstract 

We investigate the universality of correlation functions of chaotic and dis- 
ordered quantum systems as an external parameter is varied. A new, gen- 
eral scaling procedure is introduced which makes the theory invariant under 
reparametrizations. Under certain general conditions we show that this pro- 
cedure is unique. The approach is illustrated with the particular case of the 
distribution of eigenvalue curvatures. We also derive a semiclassical formula 
for the non-universal scaling factor, and give an explicit expression valid for 
arbitrary deformations of a billiard system. 
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Several criticisms were formulated many years ago concerning the applicability of random 
matrix theory (RMT) in the description of the behavior of complex quantum systems such 
as the atomic nucleus. The basic question was: should one trust the predictions of RMT at 
all if already the average density of states (Wigner's semi-circle law) does not give a good 
description of real systems? We now know that RMT describes the universal behavior of 
local fluctuation properties. These have been shown to be, in the limit of large dimensions, 
invariant for a large class of ensembles of matrices, while the average spectral density is 
ensemble dependent and therefore non-universal For a given physical system, it is 

now well established that in order to eliminate the system dependent features and observe 
universal fluctuations one should consider, instead of the original eigenenergies Ej, j = 
1,2,.. ., the unfolded spectrum 

ej = N(E 3 ) , (1) 

where 

N(E)= fp(E')dE' (2) 

is the integrated average density of states ||. By construction the new energies Sj have 
unit mean level spacing. If the system has a classical analog with n degrees of freedom 
described by the Hamiltonian H(p, q), p — (pi, . . . ,p n ), q = (qi, . . . , q n ), then to leading 
order in a semiclassical expansion the average density of states is given by the Thomas- 
Fermi approximation 

P{E) « J J d n p d n q S[E - H(p, q)} . (3) 

This expression allows for an explicit implementation of the unfolding procedure (|l|)-(|2|). 

More recently, extensions of the universal behavior of disordered as well as chaotic quan- 
tum systems have been developed to include parametric correlations and fluctuations of the 
energy levels as some real parameter X controlling the dynamics is varied (see section III-H 
of @| for a recent review). These extensions are of physical importance because many re- 
sponse functions such as the magnetic susceptibility or the conductance may be expressed as 
parametric correlations. Exactly as for the usual fluctuations computed at fixed parameter 
values, the mean properties of the flow of energy levels €j(X) when X is varied are system- 
dependent and therefore non-universal. An appropriate scaling procedure is again necessary 
in order to extract the universal behavior. 

While the average velocity of the levels is set to zero if the energy unfolding ([!]) is im- 
plemented for arbitrary X, what now characterizes the flow £j(X) is the typical slope of the 
eigenvalues with respect to X. This typical slope is not only system-dependent, but it is 
moreover parameter-dependent. In fact, the curves tj(X) as well as their average properties 
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vary when expressed as a function of a different parameter Y(X). Therefore, in order to 
have universality in the theory it is necessary to show the existence of a parametric repre- 
sentation of the elementary quantum objects satisfying two basic conditions: (i) it scales out 
system-dependent features and, (ii) it is invariant under reparametrizations. This particular 
parametric representation should moreover be unique (it should be "the" parametrization), 
since the existence of several different invariant parameters associated to X would destroy 
the universality of the correlators. 

The way to construct such a parameter is as follows. Consider the variance of the 
distribution of the parametric velocities computed at fixed X in some window (e — Ae/2, e + 
Ae/2), containing N levels, in which the statistics are evaluated 

(v 2 x ) = {l/N)Y J {de j /dXf . (4) 
j 

This function characterizes the non- universal mean properties of the flow €j (X) in the win- 
dow we are considering. More precisely, it characterizes the response of the energy levels 
to an external perturbation and, following Thouless ||, may be interpreted as a general- 
ized conductance of the system |]j|7|]. Then, to scale the spectral flow and eliminate the 
system-dependent characteristics we introduce the parameter 

V=f X 4W)dX (5) 

where Xo is some reference value. The definition of /i is universal, independent of the 
nature of X. If the system is reparametrized by introducing a new parameter Y(X), then 
\J (v x ) = J (vy) \dY/dX\. The last factor \dY/dX\ is compensated by the Jacobian of the 
transformation in Eq.(|5|) and we have 

rX , r Y 

H 



V^i) dX = sJW) dY (6) 



IX ' JY 

for any transformation Y(X). Notice that the prescription (||) yields a parametrization 
which makes the velocity variance identical to one, (v^) = 1, at all parameter values [i. 

We investigate now the uniqueness of this parametrization. It is actually possible to 
construct an infinite number of invariant parameters, according to 

ilm=M m / X (4 m ) V2m dX, 

with m — 1,2,... (all these parameters satisfy the analogue of Eq.(§)). The normalization 
constant J\f m is defined in terms of the 2m-th moment of a Gaussian distribution with 
variance one, J\f m = [(2m — l)!!] _1 / 2m . In general these parameters define different functions 
of X. However, there is one particular case for which they are all identical (and coincide 
with the simplest one m = 1 of Eq.(fj|)), and this is when the distribution of velocities is 
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Gaussian. This distribution is expected to hold for generic fully chaotic systems. It can be 
easily seen to apply to the following parametric random matrix model 



H = cosX Hi + sinX H 2 , (7) 

where H\ and H2 are two independent random matrices belonging to one of the three 
universality classes (3 = 1, 2 or 4 (orthogonal, unitary and symplectic, respectively) |§. 
Furthermore, a Gaussian distribution of velocities also holds for weakly disordered metallic 
systems, where it has been explicitly demonstrated 0. 

We have therefore established the existence of a unique, parametric-invariant scaling 
procedure for the restricted class of systems having a Gaussian distribution of velocities. 
Conversely, this distribution of velocities becomes a necessary condition for universality. In 
non-generic fully chaotic or disordered systems where the velocities are not Gaussian dis- 
tributed (like for example in strongly disordered electronic systems or banded random matrix 
models ||) the different invariant parametrizations are not equivalent, and the universality 
is lost. 

All these considerations have their analog in the usual unfolding procedure at fixed 
parameter values. The motivation to unfold the spectrum according to the prescription 
(S)"(@) - an d not according to the more "primitive" one p(Ej)Ej - is not only because it 
fixes the average mean spacing to one (both prescriptions do), but more basically because 
it makes the new spectrum ej invariant under reparametrizations of the energy. 

As an illustration let us consider two well known correlators, the velocity correlation 
function 

cM = (%M%^ + "))^, (8) 
and the distribution p{k^) of curvatures 

1 d 2 e j 

These quantities were investigated in the past for fully chaotic systems and in random matrix 
theory, as well as in disordered systems ||TU|,^|,|TT],|^,|T2|-^] . The parametric correlators were 
computed in terms of the rescaled parameter 

X = y/<&) X , (10) 

which in general is not invariant under reparametrizations and may produce non-universal 
results. For example, curvatures with respect to the parameter \x and x (k x = 
(9 2 ei /9A 2 )/(7r/3(4))) are related by 

(dej/dX) {d{vl)/dX) 
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Thus, what is expected to be universal is not the distribution of k x but the distribution of 
the particular combination given on the r.h.s. of Eq . (|TT|) . The lack of reparametrization 
invariance of k x was properly emphasized and nicely illustrated in Ref. |Tjj| . 



It follows from the definition ( JlOj ) that the parameter x coincides with fi if the function 
(v x ) is stationary with respect to X, i.e. it is independent of X. Computations done on a 
stationary spectrum having a Gaussian distribution of velocities using the parameter x are 
therefore correct in the sense that the results obtained are expected to be universal. Because 
the Hamiltonian (|7|) satisfies this property, the distribution 

p(k)=M p (i + k 2 y^i 2 (12) 



obtained from that model in Ref. jnj is the universal distribution for the curvature (here Np 
is a normalization constant). On the other hand, in the generic situation when parametric 
correlations are computed in a system where (v x ) varies with X (non-stationary spectrum), 
the use of the rescaling ([TOD leads to non-universal results and the observed distribution 



changes with X exactly as observed in Ref. [15], unless the correct parametric-invariant 



scaling (||) is used. To illustrate this point we show in Fig.l the curvature distribution of fc M 
and compare it with that of k x for the Robnik-limagon billiard (the "worst case" found in 



Ref. |y| is considered). The use of the universal parameter ji produces a dramatic change 
on the distribution. We believe that the agreement with Eq. (|T2|) will be further improved 
by going higher in the spectrum. 

The functional values of statistics involving only first derivatives with respect to the 
parameter, like in Eq.(|S]), are invariant under reparametrizations using the parameter x, and 
one would then believe that both scalings are equivalent in this particular case. However, if 
(v x ) varies with X the function looks different when plotted against x or fi, because x(X) is 
different from n(X). The use of the appropriate parametric- invariant scaling (|5|) is therefore 
necessary even for correlators involving first derivatives only. 

As for the density of states Eq.(|]), it would be desirable to have an explicit expression 
allowing for a direct computation of the non-universal function (v x ) (the generalized con- 
ductance). Such an expression may be obtained from semiclassical estimates of off-diagonal 
matrix elements by applying results of RMT || or by comparing semiclassical computations 
of the parametric density correlation function with results obtained in disordered metallic 
systems |[7l,|T7||. It is however possible to derive it from a direct semiclassical calculation 
based on the Gutzwiller trace formula. The starting point is the counting function for the 
unfolded spectrum N(e, X) = J2j @( e — e j(X)) (6 is the Heaviside function). We define the 
velocity density as 

Me) -----T 3 (^l) (13) 
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where for convenience we have replaced the delta function by an ^-smoothed Lorentzian. 
From this we obtain 



£5( e - Ci ) {de./dXf = Km2n V [/^(e)] 2 , (14) 

3 

which is our starting point for the semiclassical calculations since the average of the l.h.s. 
defines the average variance (v\). To leading order in h, N v (e,X) for a chaotic system is 
given by |T8| 

jV?( e , X) = e + % ^e iS - /h -^ /h . 
1 p -Lp 

The sum is over all the periodic orbits of the classical system, A p is an amplitude which 
depends on their stability and S p and T p are their action and period, respectively. Fur- 
thermore, t p = \dS p /de\ = \T p /p\ and the sum runs over positive and negative values of 
p. Deriving this expression with respect to X according to Eq. flT3|) , replacing in Eq. (|Tjj) , 
averaging in a small energy window and keeping only the diagonal part of the sum one 



obtains a semiclassical estimate for (v\). An analogous calculation was made in Ref. fl9| to 
compute the variance of diagonal matrix elements of an arbitrary operator A. In our case 
A = dH/dX, and we obtain 

([p " (e)]2) = h r ^>~ 2vT/m d ^ ^ 

with K D (T) = h 2 (E P \A p \ 2 5(T -T p )), and (Q 2 p ) = £ P Q 2 P \A P \ 2 / £ P \A P \ 2 with Q p = 
dS p /dX\ e . The sum runs over orbits having a period T p between T and T + dT. It can be 
shown that the vanishing of the off-diagonal terms follows from RMT, if we assume that the 
semiclassical theory can reproduce the RMT results for parametric correlations. Further- 
more, the variance of the distribution of Q p is simply proportional to the period of the orbit 
(Q 2 p) = oiT. 

The function K D (T) has also a linear dependence on the period [pi]] , K D (T) = 2 T//3. 



From this and from Eqs. flT5|) and (|TJ|) we get the final expression for the variance of the 
velocity [19| 

«> - m ■ < 16 > 

It is moreover easy to see from the previous expressions that with our assumptions there is 
no contribution to Eq. ([16]) coming from the non-universal short periodic orbits. 

The classical parameter a in Eq. (|T6"D is generally, like p, a function of X and of the energy. 
It depends on the system and the particular parametric variation under study. For example, 
for a two-dimensional billiard of area A and perimeter C, consider a general deformation 
which moves the boundary (parametrized by s, < s < C) by an amount g(s)dX in the 
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normal direction. The quantities Q v for periodic orbits with period T are given in this case 
by 

TCh 2 k 2 np 
Q P = — 2mA~^ + 2 ^E cos KM (17) 

where (/) = (l/£) J £ f(s)ds, k is the wave number corresponding to e, n p is the number 
of bounces of the orbit, and Sj and 0j denote the points of reflection and the angle of the 
trajectory with the tangent to the boundary at these points, respectively. After calculating 
(Qp) from this expression we obtain the following result for a 



a 



The constants C\ and C2 contain dynamical information about the periodic orbits, 



1 + 2E/W 



^ 2 



C2 16 



(19) 



(20) 



with 



f / \ _ (ff(s,)cosg,ff(5i +r )cosg, +T ) - (tt 2 /16)(^) 2 

/lTj (2/3)(^)-(7rVl6)( 5 ) 2 ■ 1 J 

The numerical factors in these expressions come from the ergodic phase-space averages 
(cos#) 2 = tt 2 /16 and (cos 2 fl) = 2/3. By definition, /(0) = 1 and f(r) -> for large r. The 
quantities cr 2 and nT are the variance and the average of the number of bounces of periodic 



P 



orbits with period T — > 00, respectively. In the extreme case where the correlations between 
the consecutive segments of an orbit are neglected (the uncorrelated, "random" case), then 
d = 2/3 and C 2 = 0. 

Notice from Eqs. (|T6| ) and (|I5D the energy dependence E 3 ^ 2 of the average variance of the 
velocity |16| (to leading order p = A/An for two-dimensional billiards in dimensionless units 
h = 2m = 1). This dependence holds for arbitrary deformations of chaotic two-dimensional 
billiards, but other dynamical properties of a system or other parametric variations may 
produce different energy dependences (for example, for integrable billiard systems we obtain 
(v\) oc E 2 , while for chaotic billiards with Aharonov-Bohm fluxes (v\) oc \fE JIfifl). Fig. 2 
shows the (normalized) variance of the velocity for a chaotic billiard in comparison with the 
semiclassical approximation flI§D with C\ = 2/3 and C2 = 0. A linear increase is observed, 
in agreement with the predicted energy dependence. Moreover, the slope agrees with the 
simple semiclassical estimate. The variation with energy also implies that in general it is 
not possible to simultaneously set to one the mean spacing and (v\) for arbitrary energies. 
All parametric correlators should therefore be computed in an energy window. 
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In conclusion, we have introduced a universal invariant way of scaling the parameter- 
dependent correlators in quantum mechanics. We have moreover obtained expressions which 
allow for an explicit implementation of this scaling procedure and predicted a universal en- 
ergy dependence of the variance of the velocities for deformations of a chaotic billiard. 
Contrary to the density of states which depends only on global geometrical properties of 
phase space, the derivation of Eq.([TJ]) assumes a (fully) chaotic dynamics. Moreover the 
parametric-invariant scaling procedure is unique only in the situation of generic chaotic and 
weakly disordered systems with a Gaussian distribution of the velocities. Our results, in par- 
ticular the predictions concerning two-dimensional billiards, may be tested experimentally 
in quantum dots, metallic grains or microwave cavities. 
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FIGURES 
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FIG. 1. Curvature distributions from the eigenvalues n = 2001 ton = 3000 of the odd symmetry 
class of the Robnik-limacon billiard at X = A = 0.49; dotted histogram: k x (non- invariant), 



continuous histogram: ka (invariant). Dashed curve: RMT prediction Eq.(12). 



11 



3500 




1500 1600 1700 1800 1900 2000 2100 2200 2300 2400 2500 

71 



FIG. 2. The quantity f(n) = Y7j=n {^ e j/^ r ) 2 / e j > m dimensionless units, for the eigenvalues 
of the two symmetry classes of a quarter Sinai billiard at radius r = 0.5 and no = 1500 (full lines) 
in comparison with the semiclassical approximation (dotted line). 
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